Effect of Anisotropy on the Localization in a Bifractal System 
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Bifractal is a highly anisotropic structure where planar fractals are stacked to form a 3- 
dimensional lattice. The localization lengths along fractal structure for the Anderson model defined 
on a bifractal are calculated. The critical disorder and the critical exponent of the localization 
lengths are obtained from the finite size scaling behavior. The numerical results are in a good 
agreement with previous results which have been obtained from the localization lengths along the 
perpendicular direction. This suggests that the anisotropy of the embedding lattice structure is 
irrelevant to the critical properties of the localization. 

PACS numbers: 71.30.+h, 71.23.An, 71.55.Jv 
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Recently, there has been much attention to the 

critical properties of the localization in anisotropic sys- 
tem. For a 3-dimcnsional cubic system with anisotropic 
hopping matrix, it seems now that there is a general 
agreement || that both the critical disorder (W c ) and 
the critical exponent (u) of the localization length are 
independent of direction of measurement and that such 
model belongs to the universality class of the isotropic 
Anderson model. 

On the other hand, bifractals [fflJI are constructed 
by stacking planar fractal lattices along the z-direction 
so that they are of Euclidean structure only in the z- 
direction. Anisotropy in these systems arises from the 
lattice structure itself. Therefore, it is by no means obvi- 
ous whether the critical properties obtained from the lo- 
calization lengths along the xy-plane (X xy ) are the same 
as those from the localization lengths along the z-axis 
(Az). There is a possibility that even a mobility edge, 
one of useful concepts in the localization theory, may not 
exist in this intrinsically anisotropic system. 

Therefore, in this Report, we study the critical dis- 
order and the critical exponent along fractal structure, 
i.e. along the xy-plane, for a bifractal system. Our re- 
sults are in excellent agreement with those obtained from 
A z 's, which have been reported in previous studies |Bp|, 
suggesting that the anisotropy of the embedding lattice 
structure is irrelevant to the critical properties of the lo- 
calization. In our study, one of the bifractals introduced 
in the Ref. [8] is used as the model. We consider the 
Anderson Hamiltonian given as 



(i) 



(id) 



where the random site energies are chosen from a box 
distribution of width W . The hopping energies V are 
set to 1 throughout this work and denotes the sum 
over nearest neighbor pairs of sites on the bifractal lat- 
tice. In Fig. 1(a), the lattice is schematically depicted. 
The cross section perpendicular to the z-axis is a variant 
of the Sierpinski gasket and L is the number of the frac- 
tal lattices that have been stacked. The number of the 
iteration processes for the fractal lattice is denoted by n, 



e.g. n — 2 for Fig. 1(a). This is exactly the model called 
as the Bifractal I in the Ref. [8]. The Green's function 
coupling two corner sites of the largest triangle, r and r', 
is denoted as G rjr <(n, L). Then the localization length 
along the xy-plane, X xy (L), can be defined as follows; 



X xy (L) 



= — lim 



1 



-log|G r , r <(ra,L)|, 



(2) 



where r and r' have the same z-coordinates. 

The main point of the calculation is to find G r r i (n, L)'s 
for sufficiently large value of n, which means that one 
should calculate elements of the inverse of a very large 
random matrix, (E — TL). This is essentially the same 
problem as encountered in the transfer matrix method 
for quasi- 1-dimensional systems [ [l0| ]. However a differ- 
ent recursive algorithm should be devised since we are 
considering a "quasi-2-dimensional system" . 

One can handle the problem by decimating recursively 
the amplitudes of the sites characterized by the largest it- 
eration number, as shown in Fig. 1(b). Following scheme 
is for the case of L = 1 but the extension of the method 
for L > 2 is straightforward. Let x (y) be a vector the 
elements of which are the amplitudes of the sites repre- 
sented by the solid (empty) circles in Fig. 1(b). Then a 
matrix equation for x and y can be constructed in the 
form, 



H X (E) V xy (E) \f x 
VUE) H y (E) 



(3) 



where H X (E), H y (E) and V xy (E) are matrices and con- 
tributions from remaining sites other than shown in 
Fig. 1(b) are contained in a vector z. The vector in 
the right hand side of Eq. (3) represents the fact that y 
is directly coupled only with x. Initially, when the sites 
have been indexed as in Fig. 1(b), the explicit forms of 
the three matrices are as follows; (i) H x ^j = 5ij(E — e Xi ), 



(ii) H y 
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otherwise, (iii) 



V, 
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v x 



y,36 



V and V xy _ij — otherwise. Eliminating the ampli- 
tudes of the internal sites, i.e. y, we get 



1 



V xy H y 



v xyJ 



(4) 



By performing the decimation process of Eq. (4) for every 
triangle consisting of nine sites, the number of the whole 
eigenvalue equations reduces by a factor of three. The 
3x3 matrix in the left hand side of Eq. (4) defines the 
renormalized hopping energies within the smallest trian- 
gles of the new lattice and the renormalized on-site terms. 
Since the hopping energies are modified only within the 
smallest triangle one can cast the matrix equation for the 
remaining sites again in the form of Eq. (3). It should 
be also noted that at the first step of the iteration, the 
elements of V xy are independent of E as can be seen 
from iii), while after the decimation process, i.e. Eq. (4), 
they become functions of E, in general. Therefore, the 
problem has been reduced to another eigenvalue equation 
problem on the Sierpinski gasket with the iteration num- 
ber smaller than the original by 1. Now we can iterate 
the above procedure until 3 linear equations for the am- 
plitudes of the 3 outermost sites are left in case of L = 1 . 
For general L, we have 3 x L linear equations instead 
of 3. Then the inverse of the matrix constructed by the 
coefficients of the amplitudes is the Green's function for 
the outermost sites. 

In principle, this algorithm can be iterated infinitely. 
However, one of the technical problems in real calculation 
is that some off-diagonal elements of the matrix becomes 
very small compared with, say, 0(1) diagonal elements, 
as the iteration proceeds. In general, these small ele- 
ments contain the essential information for our problem 
since they are directly related to the Green's function 
connecting two different sites far away from each other. 
For example, let's assume that the localization length 
is of order unity and the distance between two sites is, 
say, 1000 x (lattice constant). Then the Green's func- 
tion connecting the two sites is of order e" 1000 ~ 10 -430 
while the diagonal elements of the Green's function are of 
order unity. Therefore one is investigating an asymptotic 
behavior of vanishingly small matrix elements while the 
order of magnitudes of some other elements of the ma- 
trix is far much larger than them. Direct manipulation 
of such matrix leads to loss of information on the smaller 
elements. One of possible techniques to overcome this 
difficulty is to decompose the matrix into two parts as 



A = An 



(5) 



where the matrix elements of both Aq and A' are within 
the range safely handled by computers, a is a number 
and should be modified whenever the matrix is manipu- 
lated. When a matrix is decomposed as the above, the 
inverse of the matrix can be calculated by the formula 

oo 

(A + e-«A>)- 1 = V + A 1 £(- e -«A'Vr. (6) 

n=l 

For our calculation, it turns out to be sufficient to retain 
terms up to second order in e~ Q . 



We calculate A xl/ (£)'s for E = and several values of 
W in the range 4.0 < W < 9.0. In our calculation, the 
Sierpinski gasket with n = 11 (3 12 sites/cross section) is 
used and L is varied within the range 3 < L < 15. For 
a given set of parameters, configurational averages are 
performed over 4 ^ 70 different realizations to control 
the uncertainty of X xy (L) within 1 %. Periodic bound- 
ary condition is imposed in the z-direction. The results 
are shown in Fig. 2. As L increases, the renormalized lo- 
calization length, X xy (L)/ L, increases for smaller values 
of W, while for larger values it keeps decreasing. This 
implies that in the macroscopic limit, i.e., L — > oo, there 
exists a transition from an extended state to a localized 
state as W varies, e.g., from 5.0 to 6.0. Estimates of W c 
and v are obtained by fitting the data to the scaling form 



log(Ax B (L)/i) = a + b(W - W C )L X ' V . 



(7) 



where a and b are constants. Several sets of data for 
4.5 < W < 7 and 7 < L < 15, have been fitted to the 
equation (7) and finally we get W c = 5.79 ± 0.04 and 
v = 2.92 ± 0.14. The errors are the dispersions between 
different sets of data. The scaling plot with these values 
is shown in Fig. 3. Our numerical results are in good 
agreement with those obtained from the analysis of X z 's, 
i.e. W c = 5.8 and v = 3.0 ± 0.2 g for the same model. 
This supports the idea that W c and v are independent of 
the direction of measurement for this bifractal system. 

That WcS are the same along the two directions in- 
dicates that there exists a well-defined mobility edge in 
spite of the intrinsic anisotropy in this model. In addi- 
tion we can expect that the anisotropy is irrelevant to the 
localization transition in a somewhat strong sense, that 
is, the critical properties of the localization is insensitive 
not only to the anisotropy in the energy parameters but 
also to that of the embedding lattice structure. 
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FIG. 1. (a) Schematic diagram of the lattice with the it- 
eration number n — 2 and L — 5. Each vertex is a lattice 
site of the Hamiltonian given by Eq. (1). (b) Schematic dia- 
gram of the decimation process. The amplitudes of the sites 
represented by the empty circles are eliminated to modify the 
forms of the equations for the sites represented by the solid 
circles. 
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FIG. 2. Log-log plot of the renormalized localization 
length along the Sierpinski gasket as a function of L for var- 
ious values of W. The uncertainty of each data point is less 
than the symbol size. 
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FIG. 3. Scaling behavior of \og(\ xy /L) versus (W — W c ) 
L 1 /" with W c = 5.79 and v = 2.92. 
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